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Large-scale lognormal fluctuations in turbulence velocity fields 
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For several flows of laboratory turbulence, we obtain long records of velocity data. These records 
are divided into numerous segments. In each segment, we calculate the mean rate of energy dissipa- 
tion, the mean energy at each scale, and the mean total energy. Their values fluctuate significantly 
among the segments. The fluctuations are lognormal, if the segment length lies within the range 
of large scales where the velocity correlations are weak but not yet absent. Since the lognormality 
is observed regardless of the Reynolds number and the configuration for turbulence production, it 
is expected to be universal. The likely origin is some multiplicative stochastic process related to 
interactions among scales through the energy transfer. 



I. INTRODUCTION 

Turbulence spans a wide range of scales, from small 
scales in the dissipation and inertial subranges to large 
scales in the energy-containing subrange. While the small 
scales have been studied in detail, 1 the large scales have 
not thus far. There should remain unknown features that 
are crucial to understand the basic physics of turbulence. 
Here, we reveal that the rate of energy dissipation, the 
energy at each scale, and the total energy fluctuate signif- 
icantly over large scales in lognormal distributions. This 
lognormality is expected to be universal because it is ob- 
served regardless of the Reynolds number and the con- 
figuration for turbulence production. 

Landau 2 was the first to discuss such a large-scale fluc- 
tuation. The discussion was that the local rate of energy 
dissipation e should fluctuate significantly in space over 
large scales. 

Obukhov 3 also discussed the large-scale e fluctuation. 
The e fluctuation was assumed to be lognormal. This is 
because the fluctuation of any positive quantity could be 
approximated by a lognormal distribution with appropri- 
ate values for the average and standard deviation. 

The large-scale fluctuations are known to exist on some 
level, 4,5 but their details are not known. Experimentally 
or numerically, any detailed study requires long data that 
are composed of numerous ensembles. Each ensemble is 
required to cover smallest to largest scales of turbulence. 
Such long data have not been available. The situation is 
nevertheless improving for laboratory experiments, owing 
to improvements of the measurement technologies. 

We have accordingly undertaken experimental studies 
of large-scale fluctuations. 6,7 By using long velocity data 
obtained in a wind tunnel, we have shown that the large- 
scale fluctuations are always significant. 6 Their standard 
deviations are comparable to the averages. Also for a case 
of grid turbulence, we have shown that the large-scale 
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fluctuations of the energy dissipation rate, the energy at 
each scale, and the total energy are lognormal. 7 

The next problem is whether the large-scale lognormal 
fluctuations exist in other turbulent flows. For grid tur- 
bulence, boundary layers, and jets at different Reynolds 
numbers, we show that the lognormal fluctuations do ex- 
ist. The standard deviations are almost the same among 
the flows. It is thereby implied that the lognormality is 
not a mere approximation 3 but is due to some common 
process. We discuss this and other implications. 



II. EXPERIMENTS 

Table I shows experimental conditions and turbulence 
parameters for grid turbulence (Gl and G2), boundary 
layers (Bl and B2), and jets (Jl and J2). The dataset 
G2 is from our previous work. 7 We note that G2 and Bl 
as well as B2 and Jl have similar values of the microscale 
Reynolds number ReA- 



A. Wind tunnel and anemometer 

The experiments were carried out in a wind tunnel of 
the Meteorological Research Institute. We adopt coordi- 
nates x, y, and z in the streamwise, spanwise, and floor- 
normal directions. The corresponding flow velocities are 
U + u, v, and w. Here U is the average while u, v, and 
w are the fluctuations. The origin j: = i/ = z = 0mison 
the floor center at the upstream end of the test section of 
the wind tunnel. Its size was Sx = 18 m, 5y = 3 m, and 
5z = 2 m. These values of Sy and Sz remained the same 
also to x = —4 m. 

The wind tunnel had an air conditioner. If necessary, 
the conditioner was used to constrain the variation of air 
temperature. The resultant variation was ±1 °C at most 
in each experiment, where the kinematic viscosity v is 
assumed to have been constant. 

To simultaneously measure U + u and v, we used a hot- 
wire anemometer. The anemometer was composed of a 
constant temperature system and a crossed-wire probe. 
The wires were made of platinum-plated tungsten, 5 ^m 
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TABLE I: Experimental conditions and turbulence parameters for grid turbulence (Gl and G2), boundary layers (Bl and B2), 



and jets (Jl and 32). The velocity derivative was obtained 


as d x v 


[ov(x + ox) 


— 8v(x — 


Sx) — v(x 


+ lox ) + 


v(x - 25x)]/12Sx 


with 8x = U/fs- 


















Quantity 




■ 


Gl 


G2 


Bl 


B2 


Ji 


.72 


Measurement position 




m 


+1.5 


+2.0 


+ 12.5 


+ 12.5 


+ 15.5 


+ 15.5 


Measurement position 


z 


m 


1.00 


1.00 


0.35 


0.25 


0.40 


0.40 


Mean strcamwisc velocity 


u 


m s~ 


12.75 


21.16 


3.118 


11.34 


11.48 


17.38 


Sampling frequency 


fs 


kHz 


40 


70 


10 


60 


44 


70 


Kinematic viscosity 


1/ 


cm 2 s _1 


0.143 


0.142 


0.138 


0.143 


0.139 


0.139 


Mean energy dissipation rate 


(e) = 15 V ((d x vf)/2 


2 —3 

m s 


2.81 


7.98 


0.244 


12.6 


2.60 


7.52 


Kolmogorov velocity 


u K = Me}) 1/4 


m s~ 


0.0796 


0.103 


0.0428 


0.116 


0.0776 


0.101 


Rms velocity fluctuation 


(u 2 ) 1/2 


m s~ 1 


0.696 


1.10 


0.552 


2.37 


1.56 


2.34 


Rms velocity fluctuation 


(v 2 ) 1/2 


ms _1 


0.683 


1.06 


0.464 


1.96 


1.36 


2.06 


Skcwness factor 


1 3\ // 2\ 3/2 

(u )/(u ) ' 




+0.08 


+0.06 


-0.22 


-0.10 


—0.04 


— 0.04 


Skcwness factor 


(v 3 )/(v 2 ) 3/2 




—0.01 


— 0.01 


+0.01 


— 0.01 


+0.01 


+0.01 


Flatness factor 


<« 4 >/(« 2 ) 2 




3.00 


3.02 


2.69 


2.69 


2.60 


2.59 


Flatness factor 


(v 4 )/(v 2 ) 2 




2.98 


3.00 


3.06 


3.05 


3.05 


3.06 


Kolmogorov length 


n = (u 3 /(e)) 1 ' i 


cm 


0.0180 


0.0138 


0.0322 


0.0123 


0.0179 


0.0137 


Taylor microscalc 


\ = [2{v 2 )/((d x v) 2 )] 1 ' 2 


cm 


0.597 


0.548 


1.35 


0.806 


1.21 


1.08 


Correlation length 


L u — / °° (u(x + r)u{x))drj {u 2 } 


cm 


17.5 


17.9 


49.0 


43.0 


128. 


124. 


Correlation length 


L v — (v(x + r)v (x))dr/(v 2 ) 


cm 


4.46 


4.69 


6.94 


5.68 


10.2 


10.2 


Microscalc Reynolds number 


K^ = (v 2 )i' 2 \/ V 




285 


409 


454 


1103 


1183 


1603 



in diameter, 1.25 mm in sensing length, 1mm in separa- 
tion, oriented at ±45° to the streamwise direction, and 
280 °C in temperature. 



to the Gaussian values were (w 3 )/(u 2 ) 3/2 and (u 4 ) / {u 2 ) 2 . 
They tend to be sensitive to largest scale motions induced 
by either the grid, roughness, or nozzle. 



B. Flow configurations 

For the grid turbulence, we placed a grid at x = — 2 m 
across the flow passage to the test section of the wind 
tunnel. The grid was composed of two layers of uniformly 
spaced rods, with axes in the two layers at right angles. 
The cross section of the rod was 0.04 x 0.04m 2 . The 
spacing of the axes of adjacent rods was 0.20 m. We set 
the incoming flow velocity to be 12 ms" 1 (Gl) or 20 ms" 1 
(G2). The measurement position was on the tunnel axis, 
y = Om and z = 1.00 m. 

For the boundary layers, roughness blocks were placed 
over the entire floor of the test section. The block size was 
Sx = 0.06 m, Sy — 0.21m, and Sz — 0.11m. The spacing 
of the centers of adjacent blocks was Sx = Sy = 0.50 m. 
We set the incoming flow velocity to be 4ms" 1 (Bl) or 
16ms" 1 (B2). The measurement position was in the 
log-law sublayer at x = +12.5 m and y = m, where the 
boundary layer had the displacement thickness of 0.2 m 
and the 99% velocity thickness of 0.8 m. 

For the jets, we placed a contraction nozzle. Its exit 
was at x = —2m and was rectangular with the size Sy = 
2.1 m and Sz = 1.4 m. The center was on the tunnel axis. 
We set the flow velocity at the nozzle exit to be 16 ms" 1 
(Jl) or 24ms" 1 (J2). The measurement position was at 
x = +15.5 m, y = Om, and z = 0.40m. 8 

The skewness factor (tj 3 )/(u 2 ) 3 / 2 and the flatness fac- 
tor (v A ) / (v 2 ) 2 at each of the measurement positions were 
close to the Gaussian values of and 3. Here (•) denotes 
an average. These Gaussian values imply that turbulence 
had been fully developed and various eddies filled the 
space randomly and independently. 9 ' 10 Not always close 



C. Data sampling and processing 

The anemometer signal was linearized, low-pass fil- 
tered, and then digitally sampled. We set the sampling 
frequency as high as possible, on the condition that high- 
frequency noise was not seen in the power spectrum. The 
filter cutoff was at one-half of the sampling frequency. We 
obtained a long record of 4 x 10 8 data in each experiment, 
except for 1 x 10 8 data in Bl. 

The anemometer signal is proportional to the flow ve- 
locity, through the calibration coefficient that depends 
on the condition of the anemometer and thereby var- 
ied slowly in time. For individual segments of each data 
record, the length of which is fixed for the record and 
ranges from 4 x 10 6 to 2 x 10 7 data, we determined the 
coefficient so as to have the same U value. The coefficient 
within each segment is estimated to have varied by ±1% 
at most. 

We converted u(t) and v(t) at time t into u(x) and v(x) 
at position x, by using Taylor's frozen-eddy hypothesis, 
x = —Ut, which requires that (u 2 ) 1 / 2 /U is small enough. 
This requirement was satisfied in our experiments where 
(u 2 y/ 2 /U < 0.2. Even at (u 2 ) 1 / 2 /U ~ 0.3, Taylor's 
hypothesis holds valid as demonstrated by data obtained 
simultaneously with two probes separated by streamwise 
distances. 11 

Since u(t) and v(t) are stationary, u(x) and v(x) are ho- 
mogeneous, although actual turbulence in the wind tun- 
nel was not homogeneous along the x direction. Over 
small scales, fluctuations of u(x) and v(x) correspond to 
spatial fluctuations that actually existed in the wind tun- 
nel. Those over large scales do not. They are interpreted 



TABLE II: Statistics of ln5u 2 tR and ln5v 2 , R at r/L u = 0.01, 
0.03, and 0.1 as well as of In 5u R and \n5v R among segments 
with R/L u = 10 for grid turbulence (Gl and G2), boundary 
layers (Bl and B2), and jets (Jl and J2). We also show the 
number of the segments used for the analyses. 



x c , we calculate the mean energies at scale r < R as 

r x c +R/2- 



S'U 



1 f X 

1.r( x c) = -5 / 5u 2 r (x)dx, (la) 



^ r-x c +R/2-r 































Sv r (x)dx, 


Quantity 


Gl 


G2 


Bl 




B2 




Jl 




.72 




R-r 


Jx c -R/2 






Standard deviation 






















1 X 2 

ln,5 «0.01Z,u,B 


0.15 


0.14 


0.17 





15 





11 





11 














ln Su 0.03L u ,R 


0.14 


0.14 


0.16 





16 





12 





12 














ln ««S.1L U ,B 


0.16 


0.16 


0.18 





18 





14 





14 






— i- 








I 


lnu^ 

ln *«0.0H, u ,B 
ln ««S.03i u ,B 


0.31 


0.32 


0.35 





35 





25 





26 




10 2 










0.15 
0.14 


0.14 
0.14 


0.16 
0.15 






14 
15 






12 
13 






12 
13 










G2 


ln<5 «0.1L,,,B 


0.16 


0.16 


0.17 


{) 


17 





15 





15 


CM v 












lnu R 


0.22 


0.23 0.25 
Skcwncss factor 





25 





20 





20 


C\J ^ 
JO 


10 1 




v \ 






ln 6 «0.01I, u ,fl 


-0.00 


+0.06 


-0.03 


+0 


02 


-0 


03 


-0 


06 










V \ U 




lnSu l.lL u ,R 


+0.00 
-0.02 


+0.06 
+0.05 


-0.02 
-0.00 


+0 
-0 


04 
01 


-0 
-0 


05 
02 


-0 
-0 


02 
04 




10° 


/ (a) 


IOO77 

i 






In u^. 


+0.19 


+0.22 


+0.19 


+0 


08 


-0 


02 


-0 


00 












ln Sv 0.01L u ,R 
ln *«0.03£„,.R 
ln| 5«0.1Z, u ,R 


-0.01 
+0.01 
+0.00 


+0.04 
+0.05 
+0.03 


-0.13 
-0.07 
-0.07 


-0 
-0 
-0 


07 
05 
07 


-0 
-0 
-0 


03 
04 
02 


-0 
-0 
-0 


05 
05 
07 






10- 2 


10" 1 


10° 
r/Lu 


10 1 10 2 




+0.03 


+0.05 


+0.01 


+0 


01 


-0 


06 


-0 


01 















1.0 



0.5 



0.0 



Flatness factor 



lnSu 0.01L u ,R 
ln <5«0.03Z, U ,R 
ln Su 0.1L u ,R 

In u 2 r 
ln Sv, 
ln 5v, 



R 



0.03L U ,R 



In 5v Q 



ln Su 0.01L u ,R 
ln Su l.03L u ,R 
ln ««0.1L u ,B 

lntt^ 

ln Sv 0.01L u ,R 
lnSv 0.03L u ,R 
lnSv 0.1L u ,R 

N R 



2.94 
2.95 
2.97 
3.02 
2.95 
2.97 
2.99 
3.00 



3.00 
3.01 
3.01 
3.03 
3.00 
3.01 
3.01 
3.03 



3.02 
3.00 
3.02 
3.21 
3.11 
3.07 
3.03 
2.95 



3.13 
3.09 
3.09 
2.95 
3.14 
3.10 
3.01 
2.97 



3.02 
3.06 
3.02 
3.11 
3.04 
3.00 
2.97 
2.93 



3.07 
3.06 
2.98 
2.95 
2.98 
3.05 
3.03 
2.91 



Correlation coefficient against ln v R 



+0.34 


+0.35 


+0.37 


+0.37 


+0.36 


+0.36 


+0.32 


+0.33 


+0.37 


+0.35 


+0.36 


+0.36 


+0.28 


+0.29 


+0.36 


+0.32 


+0.31 


+0.33 


+0.10 


+0.12 


-0.02 


-0.05 


-0.04 


-0.03 


+0.36 


+0.37 


+0.27 


+0.30 


+0.32 


+0.33 


+0.39 


+0.39 


+0.31 


+0.33 


+0.36 


+0.36 


+0.49 


+0.49 


+0.41 


+0.41 


+0.50 


+0.49 




Number 


of segment 


3 






72600 


67600 


6360 


17560 


8120 


8000 



as fluctuations over long timescales described in terms of 
large length scales. 12 

To calculate small-scale statistics such as the mean rate 
of energy dissipation (e), we used d x v instead of usual d x u 
by assuming local isotropy, {(d x v) 2 ) = 2{{d x u) 2 ). Over 
smallest scales, the u component measured by a crossed- 
wire probe is contaminated with the w component that 
is perpendicular to the two wires of the probe. 13 The v 
component is free from such contamination. 
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III. RESULTS AND DISCUSSION 

The data records are divided into segments with length 
R. For each segment, the center of which is defined as 



FIG. 1: Statistics for grid turbulence G2 as a function of 
r/L u or R/Lu. (a) (u(x + r)u(x)) / (u 2 ) , (v(x + r)v(x)} / {v 2 } , 
and (8u 2 )/u 2 K . The arrow indicates the scale r = WOr/. 
(b) Standard deviations of ln <5«o.oil„,h hi ^ v o.03l u .r (°)j 
ln<5vo litijH (A), and lnv R (•). (c) Skewness and flatness fac- 
tors of ln6vo. 01LuiR (A), ln5«o.o3L„ i fl (°), ln^o.iL„,fl 
and lnv^ (•). (d) Skewness and flatness factors of v 2 R (o) and 
In 4 (•)• 
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with 5u r (x) = u(x+r)—u(x) and Sv r (x) — v(x+r) — v(x). 
The mean total energies are also calculated as 

I rx c +R/2 

u r(. x c) = -5 u 2 (x)dx, (2a) 

n Jx c -R/2 
^ r-x c +R/2 

v R (xc) = -5 / v 2 (x)dx. (2b) 

U Jx c -R/2 

These are coarse-grained quantities where fluctuations at 
scales smaller than R have been smoothed away. Wc 
study their statistics as a function of R. The results are 
summarized in Figs. 1-3 and Table II. 

The segment length R is given in units of the correla- 
tion length L u . This is the representative of large scales. 



Up to about 10°L„, the correlations (u(x + R)u(x)) and 
(v(x + R)v(x)) are significant [Figs. 1(a), 2(a), and 3(a)]. 
They are weak but not yet absent at about 10°L U -10 2 L U , 
which is the R range of our main interest. The scale r is 
also given in units of L u as 0.01L„, 0.03L U or 0.1L U . At 
r/L u — ► 00, 8u 2 R and Sv 2 R correspond to u 2 R and v 2 R in 
a statistical sense. 



A. Significance of large-scale fluctuations 

Figures 1(b), 2(b), and 3(b) show the standard devi- 
ations of InSv 2 R and lnv R as a function of R/L u . For 
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FIG. 2: Same as in Fig. 1 but for boundary layer B2. 



FIG. 3: Same as in Fig. 1 but for turbulent jet Jl. 



■5 



example, the standard deviation of In v R is 

<(lnt& -(lnt,*)) 2 ) 1 / 2 . (3) 

The average is over all the segments. With an increase 
of R, the standard deviations decrease, but they are not 
so small even at largest R. Thus, large-scale fluctuations 
of Sv 2 and v 2 are significant. 6,7 The same is true for 5u 2 
and u 2 (Table II). 

If we fix R/L u and r/L u , the standard deviations of 
InSulft, ln5v 2 R , and hivfj are almost the same among 
the datasets. This implies the existence of some common 
process. The exception is lnu 2 ^ that suffers from largest 
scale motions of the flow (Sec. II B). With an increase 
of r, the standard deviations of In 5u 2 R and In 6v 2 R tend 
to increase. 



B. Lognormality of large-scale fluctuations 

Figures 1(c), 2(c), and 3(c) show the skewness and 
flatness factors of ln<5v 2 R and lnvjj. For example, the 
skewness factor of In v 2 R is 

{{\nv 2 R -{\nv 2 R )f) 
((ln4-(ln4))2)3/2- ^ > 

The flatness factor of In v R is 

((ln^-qn^)) 4 ) 

({\nv 2 R - (\nv 2 R )W ( > 

At R/L u ~ 10°-10 2 , the skewness and flatness factors are 
close to the Gaussian values of and 3. Thus, large-scale 
fluctuations of 5v 2 and v 2 are lognormal. The lognor- 
mality is expected to be universal because it is observed 
in all the datasets. They are different in the Reynolds 
number or the flow configuration. Also for Su 2 R (Table 
II), we observe the lognormality. 

Figure 4 shows the probability density distributions of 
lnSv 2 R and \nv R at R/L u = 10. The distributions are 
Gaussian, at least in the observed range that extends to 
several standard deviations from the peak. 

In some datasets, deviation from the lognormality is 
seen for u 2 R (Table II). It tends toward the lognormality 
but suffers from largest scale motions of the flow (Sec. 
II B). 

The lognormal fluctuations observed here are distinct 
from those modelled by Kolmogorov. 14 While he was in- 
terested in small-scale intermittency 1 and studied 5u™ at 
small r to obtain their scaling laws, we are interested in 
large-scale fluctuations and study 5u™ R at large R. The 
scaling laws of 5u™ R are not necessary. Thus, our study 
is not necessarily concerned with the well-known prob- 
lems of Kolmogorov's lognormal model, e.g., violation of 
Novikov's inequality 15 for scaling exponents. In addi- 
tion, we are not intended to argue that the lognormality 
is exact (Sec. HIE). 



C. Correlation among large-scale fluctuations 

Table II shows the correlation coefficients of \nSu 2 R , 
lnSv 2 R , and lnu R against lnvfj. For example, the coef- 
ficient of In 5v 2 . R is 

((Infa r 2 fl -(ln^))(ln4-(ln4))) 
((ln^-(ln<5< i? ))2 ) i/2 ((ln ^_ (ln ^ ))2) i/ 2 - 1 ) 

The correlations are significant for most of the quantities. 
They tend to fluctuate synchronously. 

The datasets share interesting features. Between In u 2 R 
and lnvfj, the correlation is weak or even absent. With a 




-4-2 2 4 

(X-<X»/<(X-<X» 2 > 1/2 



FIG. 4: Probability densities P(X R ) of X R = ln5«o. iL„,-R 
(A), ln5vo. 03LutR (o), lnSvo. 1LutR (A), and lnv 2 R (•) at 
R/L u = 10 as a function of (X R - {X R ))/((X R - {X R )f) 1/2 . 
The error bars represent P(X R ) /Nj/ 2 , where Nx R is the 
number of segments that fall in the bin. The curves repre- 
sent Gaussian distributions. They are individually shifted by 
a factor 10. (a) Grid turbulence G2. (b) Boundary layer B2. 
(c) Turbulent jet Jl. 
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TABLE III: Same as in Table II but for In 5u 2 r R and In 5v* t 
at r/ij = 30, 100, and 300 as well as for In eh. 



Quantity 


Gl 


G2 


Bl 


B2 


Jl 


.72 






Standard deviation 






ln5wg R 


0.14 


0.14 


0.16 


0.16 


0.12 


0.11 


l r 2 

lnou inri „ p 

1 UUTi - ii 


0.16 


0.15 


0.17 


0.16 


0.12 


0.12 


1 A 2 

in oUgQQ^ ^ 


0.20 


0.19 


0.19 


0.17 


0.13 


0.13 


In e b. 


0.16 


0.16 


0.20 


0.19 


0.14 


0.13 




0.14 


0.14 


0.15 


0.14 


0.11 


0.11 




0.16 


0.16 


0.16 


U. lO 


0.12 


0.12 




0.20 


0.19 


0.19 


17 


0.13 


0.13 






Skewness factor 








ln Su lo V ,R 


+0.00 


+0.06 


-0.02 


+0.00 


-0.03 


-0.07 


lnSu 100r,,R 


-0.02 


+0.06 


-0.01 


+0.04 


-0.03 


-0.05 


In Su^rr* n 


— 0.02 


+0.01 


-0.01 


+0.00 


— 0.04 


-0.01 


ln e r 


-0.00 


+0.05 


-0.15 


—0.23 


-0.04 


-0.10 


lnSv 30 V ,R 


+0.01 


+0.05 


-0.09 


-0.09 


-0.04 


-0.08 




-0.00 


+0.04 


-0.06 


-0.05 


-0.03 


-0.05 


ln Sv IoO V ,R 


-0.02 


+0.01 


-0.05 


-0.07 


-0.03 


-0.06 






Flatness factor 










2.95 


3.00 


3.03 


3.13 


3.00 


3.05 


ln<5u 2 00 R 


2.97 


3.01 


3.04 


3.09 


3.02 


3.06 


1 X 2 

111 Oli'jnn,, t~j 


2.98 


2.99 


3.01 


3.07 


3.05 


3.05 


lne_R 


2.94 


2.98 


3.17 


3.35 


3.04 


2.99 


m<5t, 30,,,R 


2.97 


3.01 


3.06 


3.15 


3.02 


2.97 




2.99 


3.01 


3.03 


3.10 


3.04 


2.99 


ln<5 «300^,B 


2.97 


3.01 


3.06 


3.01 


2.97 


3.05 




Correlation coefficient against In v% 






m<5tl 30„,K 


+0.31 


+0.33 


+0.37 


+0.37 


+0.36 


+0.36 


ln5 «?00^,H 


+0.28 


+0.30 


+0.36 


+0.35 


+0.36 


+0.36 


1 A 2 

lnott 300jjR 


+0.23 


+0.25 


+0.32 


+0.33 


+0.35 


+0.35 


lnej? 


+0.35 


+0.36 


+0.23 


+0.25 


+0.29 


+0.29 




+0.39 


+0.38 


+0.29 


+0.30 


+0.32 


+0.32 




+0.49 


+0.46 


+0.36 


+0.33 


+0.33 


+0.33 


ln5 "loo^,fi 


+0.65 


+0.60 


+0.52 


+0.40 


+0.38 


+0.37 






Number of segments 








72600 


67600 


6360 


17560 


8120 


8000 



decrease of r, the correlation of ln Su 2 R against ln v R is 
increasingly strong, while that of ln Sv 2 R is increasingly 
weak. Taking also account of the observed decrease of the 
standard deviations with a decrease of r (Sec. Ill A), wc 
consider that the fluctuations originate in large r. Those 
in the u and v components are originally uncorrelated 
but become correlated in the course of their propagation 
to the smaller r. 



D. Large-scale fluctuations of small-scale quantities 

We have studied Su 2 R and Sv 2 R at r/L u = 0.01, 0.03, 
and 0.1. They are now studied at rjr\ = 30, 100, and 300. 
The Kolmogorov length rj is the representative of small 
scales. Judging from (Su 2 ) in Figs. 1(a), 2(a), and 3(a), 
the scale 30t? lies in the dissipation subrange. The scale 
100?7 lies in the inertial subrange with (Su 2 .) oc r 2 / 3 . This 
subrange extends beyond 300r? if the Reynolds number is 
high enough. We also study the energy dissipation rate 
averaged over each segment 

£r{x c ) = p / e(x)dx, (6) 

K Jx c -R/2 



with e — lbv{d x v) 2 /2. 16 It corresponds to Sv 2 at rjr\ — > 
0. 

Table III shows the results at R/L u = 10. The fluctua- 
tions are lognormal. The correlations against lnv R have 
the same features as in Table II. However, compared 
with the standard deviations at fixed r/L u in Table II, 
those at fixed r/rj in Table III tend to scatter among 
the datasets. Therefore, the fluctuations are described in 
terms of the correlation length L u that represents large 
scales, rather than in terms of the Kolmogorov length 77 
that represents small scales (see also Sec. IIIC). 

The standard deviation of ln e R is larger than that of 
ln<5v 2 0j? R , although the standard deviation of ln6v^ R de- 
creases with a decrease of r down to 30??. Probably, lne^ 
suffers from small-scale intermittency that causes a broad 
e distribution. 1 

For e R in the dataset Bl and B2, we see deviation 
from the lognormality. This is achieved at the larger R. 
In fact, at R/L u = 100, the skewness factor is +0.03 (Bl) 
or +0.00 (B2). The flatness factor is 3.04 (Bl) or 3.00 
(B2). 

E. Origin of large-scale lognormal fluctuations 

Regardless of the Reynolds number and the flow con- 
figuration, the energies <5w 2 and 5v% at scale r, the total 
energies u 2 and v 2 , and also the energy dissipation rate 
£ fluctuate over large scales in lognormal distributions. 
Their standard deviations are almost the same among 
the flows (Sees. Ill A, IIIB, and HID). Hence, some 
common process is at work, which is discussed here. 

The observed fluctuations originate in turbulence it- 
self. They are too strong to be an experimental artifact, 
i.e., some variation in the condition of the wind tunnel 
or the measurement system. In addition, if such a vari- 
ation were responsible for the fluctuations, they should 
exhibit a correlation between u 2 R and v R that is actually 
not significant at all (Sec. IIIC). 

A lognormal fluctuation is induced by some multiplica- 
tive stochastic process, e.g., a product of many inde- 
pendent stochastic variables, Qn — Yln=i 9™ ■ ^ ts l°S a " 
rithm is a sum of many independent stochastic variables, 
ItlQn — ^2n=i m <?ni which is Gaussian according to the 
central limit theorem. For the lognormal fluctuations ob- 
served here, the process is related to interactions among 
scales through the energy transfer. 7 While the mean en- 
ergy transfer is to a smaller scale and is significant be- 
tween scales in the inertial subrange alone, the local en- 
ergy transfer is either to a smaller or larger scale and is 
significant between all scales. 6 ' 17-19 Thus, any scale in- 
teracts with itself and with many other scales. 

The fluctuations observed for different scales r tend to 
be synchronous (Sec. IIIC). This is consistent with the 
above explanation that the fluctuations are induced by 
interactions among scales. 

The rate of the energy transfer does not have a lognor- 
mal fluctuation because the rate could change its sign. 
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Although the fluctuation of this rate is important, it is 
beyond the reach of our experimental study. 

To apply the central limit theorem, there exist various 
forms of sufficient but not necessary conditions. 20 The 
essence is that the stochastic variables are all alike and no 
few of them dominate the others. It is also required that, 
among the variables, correlations are absent exactly or at 
least practically. These requirements have to be satisfied 
by the interactions among scales that are to reproduce 
the observed lognormality. 

The interactions are also required to involve a wide 
range of scales because the lognormal fluctuations are ob- 
served up to large R and down to small r. Large scales 
are especially important as implied by the observed fea- 
tures of the standard deviations and correlations (Sees. 
Ill A and niC). 

The central limit theorem does not determine the dis- 
tribution at its tail, which might deviate from the lognor- 
mal distribution and even might depend on the Reynolds 
number or the flow configuration. Whether such a tail is 
prominent depends on the case. An extreme case could 
be the observed deviation from the lognormality of u 2 R in 
some datasets (Sec. IIIB). Nevertheless, the lognormal- 
ity is important at least as an idealized model because it 
exists only if the logarithm is subject to the central limit 
theorem. The lognormality also accounts for most of the 
distributions in the observed range (Fig. 4). To observe 
the distributions over the wider range, much longer data 
are necessary. 



F. Scale range of lognormal fluctuations 

The lognormal fluctuations exist at R/L u ~ 10°-10 2 , 
where the correlations (u(x + R)u(x)) and (v(x + R)v(x)) 
are weak but not yet absent [Figs. 1(a), 2(a), and 3(a)]. 
Although the correlations there exhibit different profiles 
among the flow configurations, this difference is not rel- 
evant to the existence of the lognormality. Outside the 
R range, the correlations are either too strong or absent. 
The lognormal fluctuations do not exist. 

At R/L u < 10°, the skewness and flatness factors of 
In Sv 2 R and In v R deviate from the Gaussian values [Figs. 
1(c), 2(c), and 3(c)]. The fluctuations are not lognormal 
and no longer explained by a multiple of independent pro- 
cesses. This is because the correlations (u(x + R)u{x)) 
and (v(x + R)v(x)) are too strong. The observed devia- 
tion from the lognormality is different among the dataset. 
Note that, among the figures, different are the ordinates 
for the skewness and flatness factors. Since this R range 
for small scales has been studied in detail, 1 we do not 
study it any more. 

At R/L u > 10 2 , the correlations (u(x + R)u(x)) and 
(v(x+R)v(x)) are absent. Then, any segment with length 
R is divisible into many subsegments that are indepen- 
dent of one another. The length i?* of these subsegments 
is such that L u -C R* <C R- By using a sum over the sub- 



segments, we write sr, 5u 2 r , Sv 2 R , u R , and v R as, e.g., 
R R/R, 

= ^ E (7) 

n— 1 

Here x Cn is the center of the n-th subsegment. The sum 
is subject to the central limit theorem and thereby fluc- 
tuates in a Gaussian distribution. 6 ' 21 Figures 1(d), 2(d), 
and 3(d) show the skewness and flatness factors of v R , 
which in fact approach to the Gaussian values with an 
increase of R/L u . 22 ' 23 In other words, when R/L u is large 
enough, quantities such as Rv R are additive. 6 The value 
of an additive quantity for a region is the sum of its values 
for the subregions that are independent of one another. 24 
Such additive quantities are well known in thermodynam- 
ics and statistical mechanics, for which correlations are 
absent over scales of interest. 



IV. CONCLUDING REMARKS 

For several flows of laboratory turbulence, we obtained 
long records of velocity data. These records were divided 
into numerous segments with length R. In each segment, 
we calculated the mean rate of energy dissipation sr, the 
mean energies 5u 2 R and 5v 2 R at scale r, and the mean 
total energies u 2 R and v R . Their values fluctuate signif- 
icantly and synchronously. The fluctuations are lognor- 
mal. Since the lognormality is observed regardless of the 
Reynolds number and the configuration for turbulence 
production, it is expected to be universal. 

The lognormal fluctuations are of large scales, R/L u ~ 
10°-10 2 , where the correlations (u(x+R)u(x)) and {v(x+ 
R)v(x)) are weak but not yet absent. This R range had 
not been studied, although it lies in between the ranges 
that had been studied in detail. 

We explained the lognormal fluctuations by applying 
the central limit theorem to the logarithm of a multiplica- 
tive stochastic process, which was related to interactions 
among many scales through the energy transfer. To ap- 
ply the central limit theorem, it is required that the inter- 
actions are all alike and do not have mutual correlations. 
In addition, features of the observed fluctuations require 
that the interactions involve a wide range of scales. Large 
scales are especially important. 

The first discussion on large-scale fluctuations was by 
Landau, 2 who argued against the universality of small- 
scale statistics discussed by Kolmogorov. 25 Kolmogorov 
discussed that small-scale statistics depend on v and (e) 
alone and are thus universal. Landau discussed that 
small-scale statistics are not universal because the large- 
scale e fluctuation is significant and is not universal. In 
practice, large-scale fluctuations are universally lognor- 
mal. Their standard deviations are also almost univer- 
sal. Even if large scales affect small scales through the e 
fluctuation, 3,26 it does not necessarily violate the univer- 
sality of small-scale statistics. However, past experiments 
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suggested that some small-scale statistics are not uni- 
versal and depend on the large-scale flow. 6 ' 11, 27-29 Since 
the Reynolds numbers of these studies might not have 
been high enough to distinguish between small and large 
scales, 30 ' 31 it is desirable to study the effect of large scales 
on small scales at the higher Reynolds numbers. 

Finally, we note that the large-scale lognormal fluctu- 
ations are not restricted to those of turbulence. They 
should also exist in other nonlinear systems that have 
many degrees of freedom. This is because, regardless of 
the process details, lognormality exists only if the loga- 
rithm is subject to the central limit theorem. However, so 
far as we could see, no example is known. It is desirable 



to explore the lognormal fluctuations in various systems, 
promisingly over length scales or timescales where corre- 
lations are weak but not yet absent. We expect that the 
large-scale lognormal fluctuations turn out to be ubiqui- 
tous, far beyond those of turbulence observed here. 
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